Point defects on graphene on metals 
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Understanding the coupling of graphene with its local environment is critical to be able to integrate 
it in tomorrow's electronic devices. Here we show how the presence of a metallic substrate affects the 
properties of an atomically tailored graphene layer. We have deliberately introduced single carbon 
vacancies on a graphene monolayer grown on a Pt(lll) surface and investigated its impact in the 
electronic, structural and magnetic properties of the graphene layer. Our low temperature scanning 
tunneling microscopy studies, complemented by density functional theory, show the existence of a 
broad electronic resonance above the Fermi energy associated with the vacancies. Vacancy sites 
become reactive leading to an increase of the coupling between the graphene layer and the metal 
substrate at these points; this gives rise to a rapid decay of the localized state and the quenching of 
the magnetic moment associated with carbon vacancies in free-standing graphene layers. 

PACS numbers: 73.22.Pr, 73.20.Hb, 68.37.Ef, 71.15.Mb 



Exciting properties were supposed for graphene since 
long time ago However, it was not till 2004 [2] that 
graphene ceased being a theoretical chimera to become 
the object of desire of the scientific community. In just 
few years, most of these extraordinary properties have 
already been demonstrated [3-5] and many others are 
emerging as a result of the tremendous experimental and 
theoretical efforts devoted to this material [sl, 0] • As a 
consequence, graphene has undoubtedly become one of 
the most promising candidates to play a key role in future 
technology. Many of the experimental efforts have been 
invested in growing larger and higher quality graphene 
layers and also in understanding and controlling the cou- 
pling of graphene with other materials, in particular with 
metals, a must to incorporate graphene to real devices. 
Epitaxial graphene on metals represents an ideal route to 
fulfill both requirements. Highly perfect graphene sheets 
can be grown on various metals [8-10] and very recently 
macroscopic-sized graphene films have been grown and 
subsequently transferred to arbitrary substrates . 
The interaction of graphene layers with the different 
metallic substrates strongly depends on the metal itself 
13|, [l^ , but according to its strength two main groups 
can be identified: strongly and weakly interacting sys- 
tems [13]. While properties of graphene monolayers be- 
longing to the first group can be quite different with re- 
spect to the free-standing case [15|, in the weakly inter- 
acting systems the electronic structure of ideal graphene 
is basically preserved, as revealed by the experimental 
observation of Dirac cones comparable to those of per- 
fect graphene [l6|, [l^ . 

Till date, graphene on metals research has basically 
focused on pristine graphene monolayers adsorbed on 
metal surfaces, which are at present quite well under- 
stood [IGI, [iTj . In contrast, despite the pressing chal- 
lenge to tailor graphene layers in order to fully exploit 



graphene's potential, very little is known about the influ- 
ence of the metallic substrate in the properties of modi- 
fied graphene layers [19]. Here, we show that while prop- 
erties of pristine graphene adsorbed on a weakly inter- 
acting metal like Pt(lll) are reasonably preserved, the 
situation is dramatically different when point defects are 
introduced in the graphene layer. 

Our starting point is a perfectly clean graphene mono- 
layer adsorbed on Pt(lll), one of the weakest interact- 
ing graphene- metal systems [lol, [l3 • As a result of such 
a weak interaction, the graphene layer presents several 
orientations with respect to the Pt(lll) surface, giv- 
ing rise to various moire patterns 0, [Til, [20|, |21|. In 
addition, recent photoemission experiments have shown 
the existence of tt bands presenting linear dispersion, 
the same Fermi velocity as in free standing graphene 
(FSG) and a Dirac point slightly shifted to ~ +300 meV 
[rz]. We formed a complete graphene monolayer on the 
clean Pt(lll) surface by chemical vapor deposition of 
ethylene in ultrahigh vacuum (UHV) environments at 
temperatures above 1275 K. Our experiments were per- 
formed at 6 K using a home-made low temperature scan- 
ning tunneling microscope (LT-STM) operating in UHV 
[22]. Figure la shows a large STM image of an atomi- 
cally perfect graphene/Pt(lll) surface, where two moire 
patterns can be identified. In the upper-left corner an 
(a/2T X a/2T)^11° moire is found, while the rest of the 
image shows a 3 x 3 moire (both periodicities with respect 
to the graphene lattice). Such a 3x3 superperiodicity is 
due to a 19.1° angle between graphene and the Pt(lll) 
surface, corresponding with a lattice misfit of 0.60% (2Qj . 

All the experimental results shown here are basically 
independent of the moire periodicity, thus, for the sake 
of simplicity and to facilitate a straightforward compar- 
ison with theory, we will restrict ourselves to the 3x3 
moire (Fig. lb), one of the most frequently found in the 
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FIG. 1. (color online) a) STM image of the pristine 
graphene/Pt(lll) surface showing two different moire struc- 
tures, b) Zoom in of the 3x3 region highlighted in a). Sample 
bias: 50mV, tunneling current: l.OnA for a) and b). STM 
data were measured and analyzed with WSXM [42.] . c) Simu- 
lated STM image of 3x3 graphene/Pt at V=+100 mV. d) 
Calculated band structure of the 3x3 Graphene/Pt moire 
(black lines) and pure graphene (red lines). Pure graphene 
bands have been shifted by +410 meV. e) STS measurement 
of the LDOS of pristine graphene on Pt(lll).f) Theoretical 
DOS of the 3x3 moire projected on the p states of C atoms 
(s and d contributions are negligible in this energy range). 

graphene/Pt (111) system pjQ. Information about the lo- 
cal density of states (LDOS) of the sample was obtained 
with atomic precision by measuring at 6 K differential 
conductance (dl/dV) spectra in open feedback loop mode 
using the lock-in technique with frequency 2.3 kHz and ac 
modulation of 1 mV. Our dl/dV spectra measured in the 
pristine graphene surface show a clear dip at ~ +300 mV 
accompanied by a Fshaped rise at both sides, see Fig. le. 
Recent scanning tunneling spectroscopy (STS) experi- 
ments measured on a partially covered graphene/Pt (111) 
surface, have reported dl/dV spectra with a slight dip at 
+150 mV [23], which was related to the unoccupied sur- 
face state existing in Pt(lll) [2^. However, the clear V 
shape of our spectra around the dip and its location at 
+300 mV, coinciding with the Dirac point energy the- 
oretically predicte d [ 13j and estimated from photoemis- 
sion experiments [17|, make us believe that the dip we 
observe at +300 mV is associated with the position of 
the Dirac point in the graphene/Pt (111) surface. An- 
other piece of evidence comes from the results of our 
density functional theory (DFT) calculations [25] based 
on VASP [26]. We have studied the 3x3 moire using 
an empirically corrected version of the PBE functional 
[27] that includes the effects of van der Waals interac- 
tions [28] . The calculated bands of the system are shown 
in Fig. Id (in black). In the same figure we also in- 
clude the bands of pure graphene (in red) shifted by 
+410 meV to make them overlap with their equivalent in 
the moire. Consistently, the calculated DOS of Fig. If, 
which can be directly compared to the experiment in Fig. 
le, presents a V-shaped minimum at that same energy. 



Our theoretical results capture the structural properties 
of the graphene/Pt (111) system very well, in opposition 
to other methodologies that include disp ersion interac- 
tions into standard DFT functionals [29]. In particular, 
our calculated average distance between graphene and 
the uppermost Pt layer is 3.35 A, whereas LEED/LEEM 
experiments yield 3.30 A[17]. 

Experimental STM images of graphene adsorbed on 
Pt(lll) (Fig. lb) can be understood in more depth 
with the help of theory. In our approach, we use a non- 
equilibrium Green's function formalism to evaluate the 
currents [30|, using the OpenMX code [31] to map the 
hamiltonian into a local orbital basis, and an idealized 
Pt apex with a single dz^ orbital to represent the mi- 
croscope tip. This model produces atomically resolved 
images (see Fig. Ic) in good agreement with the exper- 
imental results in Fig. lb. Surprisingly, the brightest 
features in the calculated image correspond to C atoms 
that are lowest in the graphene sheet. Notice that in 
the 3x3 moire the lattice mismatch is very small, and 
the topographic corrugation, measured as the height dif- 
ference between the lowest and highest C atoms, is only 
0.02 A. In this case, the observed anticorrelation between 
the simulated image and the atomic topography implies 
that the STM corrugation is not a geometrical effect but 
instead it can be explained as a purely electronic effect. 

To atomically tailor the graphene layer adsorbed on 
Pt(lll), we have generated point defects on it by irradi- 
ating the surface with 140 eV Ar+ ions, which are known 
to mainly produce single C vacancies on a graphite sur- 
face [32]. Our LT-STM images show that the previously 
pristine graphene sheet presents now a number of almost 
identical bright features, associated with the amount of 
Ar+ ions reaching the surface, see Fig. 2. The exact 
shape of these features slightly depends on its location 
with respect to the moire pattern, as shown in Figs. 2b- 
e. We see an elongated protrusion occupying an exten- 
sion of two honeycomb lattice parameters, surrounded 
by a non trivial pattern of high LDOS intensity extend- 
ing less than Inm. In ideal graphene, the presence of 
a point defect should generate short-wavelength modu- 
lations of the LDOS with (^3 x V^) R30° (R3 in the 
following) periodicity due to intervalley scattering [33]. 
In well-decoupled graphene layers, as epitaxial graphene 
on SiC or HOPG surfaces, such R3 patterns associated 
with atomic-size impurities have been observed by means 
of low bias STM images, which are a measure of the 
LDOS at E^, [3i,[3ii3- 

2D Fourier transforms of our 
low bias STM images measured in samples with a high 
enough density of generated atomic vacancies, also al- 
lowed us to measure such R3 modulations and to esti- 
mate a Fermi wave- vector kp ^ 0.5nm~^, which is con- 
sistent with Ei:)=+300meV and a Fermi velocity lO^m/s 
as recently measured by photoemission [13, [111 . 

Contrary to the case of the graphite surface [lH, 
we cannot directly identify point defects generated by 
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FIG. 2. (color online) a) STM topography, measured at 6 K, 
showing the graphene/Pt(lll) surface after the Ar^ irradi- 
ation. Sample bias: -30mV, tunneling current: 0.8nA. b-e) 
Zoom-ins showing the small variation in the bright features 
shape for different moire positions. All bright features on the 
image can basically be identified with one of these 4 images. 



Ar sputtering in the graphene/Pt(lll) surface with the 
three fold patterns predicted for single C vacancies on a 
FSG layer Thus, we have to bring into play DFT 

calculations in order to unravel the nature of the atomic 
point defects observed in this system. In our simulations, 
we use a supercell with 2x2 units of the 3x3 moire and 
a single C atom removed j25|. For this structure, several 
adsorption sites are possible, depending on the relative 
location of the vacancy and the underlying Pt surface. 
We have studied some of these possibilities, finding that 
the best match to STM experiments occurs for structures 
which reconstruct similarly to single vacancies on isolated 
graphene sheets [ssj . In the reconstruction of free stand- 
ing graphene, two of the 3 under coordinated C atoms 
surrounding the vacancy move closer to each other and 
become weakly bonded, forming a pentagon ring, and the 
third undercoordinated atom moves out of the graphene 
plane by only ~ O.lA. In graphene/Pt(lll), the third 
atom and one of its neigbours move out of the plane and 
towards the Pt by ~ 1 A, and form two new chemical 
bonds with Pt atoms of the surface (Figs. 3a and b). 
Besides, the average distance between the graphene sheet 
and the topmost Pt layer has decreased by ~ 0.08 A. It 
is thus clear that the previously inert graphene layer has 
become very reactive. The strong interaction between 
the graphene with vacancy and the metal also results in 
a quench of the possible magnetic moment of the system, 
with the relaxed structure being non-magnetic. The sim- 
ulated STM images present a pattern like that in Fig. 3c. 
The brightest features are a heart-shaped, elongated pro- 
trusion and a small, oval protrusion next to it. Another 
characteristic feature is a small dark area right next to 
the elongated protrusion. All these structures can be 
recognized in the experimental images (Fig. 3d). By 
comparing the simulated image with the position of the 
atoms (Fig. 3c) we can associate the heart-shaped fea- 
ture to the region where the two weakly bonded C atoms 




FIG. 3. (color online) Relaxed single carbon vacancy in a 6x6 
unit cell (4 cells of the 3x3 graphene/Pt(lll) moire ), with 
the 2 C atoms that move towards the Pt highlighted: a) top 
view, b) side view, c) calculated constant current STM image 
at V = +100 mV with the ball and stick atomic atomic model 
superimposed, d) Experimental STM image of a vacancy from 
panel 2a. 



and the vacancy itself are located, the dark region to the 
out-of-plane C atoms, and the small, oval feature to the 
other atom in the pentagon ring which belongs to the 
same sublattice as the weakly bonded pair. The modu- 
lation of the charge density revealed by the STM image 
decays fast. However, we do not try to compare it quan- 
titatively with the experiment, since our unit cell is still 
small to model an isolated vacancy, and we believe the 
intensity decay can be affected by interference between 
neigbouring cells. 

Introducing single C vacancies in an isolated graphene 
layer has a profound impact in its properties. Accord- 
ing to many theoretical studies 39|-41|, they give rise to 
quasi-localized states at the Fermi level, which can be as- 
sociated with the generation of local magnetic moments 
around the C vacancies and produce a strong reduction 
of charge carriers' mobility. Such theoretical expecta- 
tions were recently confirmed by some of us by STS ex- 
periments on C vacancies on the graphite surface, which 
revealed the presence of a very sharp resonance at the 
Fermi energy extending more than 3 nm away from each 
single C vacancy [32|. It is thus crucial to understand 
how the coupling with the metallic substrate will affect 
the properties of such C vacancies. 

Our STS measurements on C vacancies in 
graphene/Pt(lll) show a strong increase of the 
LDOS starting at the Dirac point and reaching a 
maximum around +500 meV. As it can be seen in 
Fig. 4a, this is reflected in our dl/dV spectra as a 
broad electronic resonance centered at +500 meV with 
a FWHM ~150 meV. dl/dV spectra acquired on C 
vacancies located in different positions of the 3x3 moire, 
show that all C vacancies present this broad electronic 
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FIG. 4. (color online) a) STS measurements of the LDOS on 
a C vacancy (blue) and on pristine grapliene/Pt(lll) (red). 
dl/dV were taken consecutively at 6 K with the same micro- 
scopic tip. b) dl/dV spectra measured on different C vacan- 
cies, numbered 1 to 5. c) dl/dV map, measured at a sample 
bias of +500mV. 1 to 5 point to same vacancies as panel b). 



resonance at +500 meV independently of its position 
inside the moire, which points to a small influence of 
the moire superstructure on its properties. The only 
influence detected was the variation of the resonance 
height for C vacancies on different moire positions, see 
Figs. 4b-c. Analogous results were also obtained for 
different moire patterns. Since the electronic resonance 
associated with any vacancy lays at a fixed energy of 
+500 meV, it is possible to map its spatial extension 
by conductance images at +500 mV. Figure 4c shows 
one of such conductance images measured on a 15x15 
nm^ region where 11 single C vacancies have been 
generated by Ar+ irradiation. As it can be observed, the 
localized state only extends 0.6 nm away from the center 
of the vacancies. The conductance map also reflects 
the resonance height variation found for vacancies in 
different moire locations. Our calculations show that 
this peak is mainly caused by the 7r-electrons of the two 
weakly bound atoms in the graphene reconstruction. 

Thus, our results show that the electronic properties of 
C vacancies in a graphene layer on top of a metal strongly 
differ from the ones found for C vacancies in a well decou- 
pled graphene layer as it is the graphite surface |32j. As 
shown by our DFT calculations these differences can be 
explained as due to the increase of the interaction with 
the metallic surface after the C vacancy formation. 

In summary, our findings demonstrate that even in 
weakly coupled graphene/metal systems the presence of 
the metal has to be seriously taken into account in order 
to controllably tune graphene properties by locally mod- 
ifying its structure. We have shown that the interaction 
with the metal strongly increases when single C vacan- 
cies are introduced in the graphene layer, discarding the 
appealing possibility of inducing magnetic properties on 
this graphene system by removing C atoms. 
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